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Lecture  1 


PRE-ASYMPTOTIC  COMBUSTION  REVISITED 

The  description  of  reacting  systems  can  be  simplified  when  the  so-called 
activation  energy  is  large;  the  notion  is  an  old  one,  but  its  full  power  is 
only  released  by  modern  singular  perturbation  theory.  More  than  forty  years 
ago,  Frank-Kamenetskii  introduced  approximations  based  on  large  activation 
energy  to  construct  a  thermal  theory  of  spontaneous  combustion,  and  we  shall 
start  there.  His  problem,  which  neglects  the  fluid-mechanical  effects  of 
main  concern  to  us ,  focuses  attention  on  the  reaction  and  thereby  acts  as  a 
precursor  for  the  lectures  that  follow.  The  problem  and  its  generalizations 
have  been  the  happy  hunting  grounds  of  mathematical  analysts  for  many  years, 
but  it  was  not  until  quite  recently  that  a  complete  description  of  the 
ignition  and  explosion  processes  was  made  available  by  Kapila  and  Kassoy 
(working  separately)  through  activation-energy  asymptotics,  the  main  theme 
of  these  lectures. 

1.  Ignition 

Let  us  suppose  that  a  combustion  system  has  a  characteristic  temperature 

,Tc  and  that  the  heat  generated  by  reaction  can  be  expressed  as  a  function  of 

T  in  the  Arrhenius  form 
c 

-9/T 

qe  C.  (1) 

This  function  has  an  inflection  point  at  T  *  0/2  and  its  second  derivative 

c 

is  positive  for  smaller  values,  where  the  graph  is  accordingly  concave  upwards 
(figure  1).  Suppose  also  that  the  heat  loss  by  conduction  and  convection  has 
the  linear  form 


k(Tc-Tf) 

where  is  the  ambient  temperature. 


(2) 


-1.2- 


The  system  can  only  be  in  equilibrium  if  the  heat  generated  (l)  is  equal 
to  the  heat  lost  (2).  The  parameters  q  and  k  do  not  then  play  independent 
roles  but,  rather,  it  is 

V  =  q/k  (3) 

that  is  relevant.  This  ratio,  which  will  be  called  the  Damkohler  number  (cf. 
later  lectures),  can  be  altered  by  changing  the  parameters  of  the  system. 

It  is  apparent  from  figure  1  that,  for  9  >  UTf,  there  is  either  one 
solution  or  there  are  three  solutions,  depending  on  the  value  of  V.  If  V 
is  increased,  the  straight  line  rotates  about  its  end  point  in  a  clockwise 
direction’,  and  we  can  identify  the  transitions  from  1  solution  (cold)  to  3 
solutions  to  1  solution  (hot).  It  is  this  second  transition  that  is  our 
concern  in  this  lecture.  The  state  of  the  system,  represented  by  a  cold  point 
such  as  C,  moves  towards  I  as  the  Damkdhler  number  is  increased,  and  then  must 
Jump  to  a  hot  point  such  as  tH  at  the  transition.  This  jump  is  called 
ignition. 

Ignition  is  ubiquitous  in  combustion  systems;  it  can  generally  be 
attributed  to  the  nonlinear  dependence  of  heat  generation  on  temperature  and 
the  essentially  linear  dependence  of  heat  loss.  The  precise  nature  of  the 
phenomenon  can  only  be  determined  by  detailed  analysis,  though  the  results 
of  different  calculations  carried  out  by  activation-energy  asymptotics  often 
bear  a  strong  family  resemblance.  They  are  characterized  by  the  following 
elementary  example,  introduced  by  Frank -Kamenetskii. 

2 .  Spontaneous  Combust ion 

Consider  the  boundary -value  problem 

d2T/dx2  -  -  Pe"9/T  for  |x|  <  1,  T  «  Tf  at  x  «  tl.  (U) 

Heat  conduction  in  the  infinite  slab  is  balanced  by  heat  generation  due  to  the 


reaction.  Depletion  of  the  reactant  has  been  ignored,  so  that  the  reaction 
rate  depends  only  on  temperature,  as  in  section  1. 

Such  models  have  been  used  for  many  years  to  explain  spontaneous 
combustion,  the  auto-ignition  that  occurs,  for'  example,  in  large  volumes 
of  damp  organic  material.  Nondimens  ionalizat  ion  of  the  heat -conduct  ion 
Laplacian  makes  the  Damkohler  number  V  proportional  to  a  ,  where  a  i3 
a  length  characteristic  of  the  volume  (here  the  semi -thickness  of  the  slab). 
Thus,  an  increase  in  V  may  be  achieved  by  increasing  the  volume  and,  a3  we 
shall  see,  this  can  lead  to  ignition. 

We  seek  a  solution  of  the  problem  (h)  that,  as  3  -*■  ■»,  deviates  only 
by  0(9-^)  from  the  uniform  state,  i.e. 

T  =  Tf  +  0-1T^  +...  with  <j>  *  ( 1/T ) x •  (5) 


This  leads  to 


?  ?  *  -9/X  o 

d  <fr/dx,  =  -5e9  with  5  *  P9e  1  /T^, 


(6) 


an  equation  first  obtained  by  Frank -Kamenetskii,  and  the  boundary  conditions 


$  =  0  at  x  *  tl. 


(T) 


Here  5,  the  scaled  Damkohler  number, is  assumed  to  be  0(l). 

The  perturbation  <P  achieves  its  maximum  ($m)  at  the  midpoint  x  *  0 
and  so  may  be  written 


♦  / 2 


$  *  2  in[e  m  sech(cx)]  with  cc  *  j<5e 


m 


(8) 


The  boundary  conditions  then  imply  that 


t——  .  ♦  /« 

/6/2  ■  e  cosh  ^(e  m  2), 


(9) 


which  defines  the  maximum  temperature  in  terms  of  the  parameter  5.  The  result 


is  shovn  in  f ignore  2,  displaying  the  phenomenon  of  ignition.  For  <$  less  than 
the  critical  value 


5c  =  0.878,  (10) 

there  is  a  steady  solution,  in  fact  there  are  two  solutions;  whereas  for 
6  >  there  is  no  solution  of  the  type  (5),  and  in  fact  none  at  all.  The 
absence  of  a  steady  state  for  supercritical  values  of  5  implies  that,  with 
unsteady  effects  included,  the  temperature  will  increase  without  hound  when, 
for  example,  the  system  is  initially  in  a  uniform  state  T  *  T^.  In  practice, 
the  increase  is  limited  by  depletion  of  the  reactant,  an  effect  that  is 
ignored  here.  Of  the  two  solution  branches  for  5  <  6^,  the  upper  one  is 
believed  to  be  unstable,  though  this  has  never  been  proved. 

It  is  a  general  characteristic  of  ignition  that  it  is  associated  with 
0(9-1)  perturbations  of  the  frozen  solution,  i.e.  the  solution  obtained  for 
0  =  0.  This  is  certainly  true  for  the  diffusion  flames  treated  in  lectures 
8  and  9  (see  section  8.6);  the  details  differ  from  those  presented  here, 
but  the  essential  ideas  do  not. 

3.  Homogeneous  Explosion 

So  far  we  have  inferred  ignition  from  a  steady-state  theory.  The 
phenomenon  itself  is  inherently  unsteady,  and  certain  aspects  of  the  unsteadi¬ 
ness  deserve  examination.  To  that  end,  it  is  useful  to  consider  first  the 
spatially  homogeneous  initial  value  problem 

dT/dt  »Pe"0/T,  T  »  Tf  for  t  »  0.  (11) 


There  is  no  value  of  V  for  which  a  steady  state  can  be  attained;  the 
problem  is  always  supercritical.  The  physical  reason  is  that  no  heat- 
loss  mechanism,  such  as  conduction  to  the  boundaries  (section  2),  exists. 


-lo¬ 


in  terns  of  the  exponential  integral 

y  u 

Ei(y )  =  f  —  du. 


(12) 


this  has  the  exact  solution 


Vt  -  9[f(T  )-f(T)]  with  f(y) 


for  any  value  of  0.  Since  the  function 


f(y)  *  (y/0)2e9/y  +...  as 


Ei(0/y)-(y/0)e0/y  ' 
f  has  the  symptotic  expansion 

0  -*■ 


(13) 


(1U) 


it  follows  that 


.  .  2  0/T-6/T 

t  *  te~(T/Tf)  e  I+...  with  t  =  5t; 


(15) 


here 

£  *  1.  (16) 

s 

m 

Deviations  from  the  initial ‘state  of  the  form  (5)  are  therefore  described  by 


$  *  -in(tg-t) ,  (17) 

so  that  t  is  the  time  to  explosion  (i.e.  the  time  that  T  takes  to  deviate 
from  its  initial  value  T^.  by  more  than  0(0  ^)).  The  behavior  of  T  in 
some  small  neighborhood  of  t  is  called  thermal  runaway. 

Thermal  runaway  terminates  what  is  known  as  the  induction  phase,  the 
problem  with  which  pre-asymptotic  theory  was  almost  exclusively  concerned. 

To  go  substantially  further,  modern  asymptotic  methods  must  be  used,  as  will 
be  discussed  in  these  lectures.  To  analyze  the  so-called  explosion  phase 
that  follows  induction,  we  return  to  the  expansion  (15),  and  introduce  a 
fast  time  t  given  by 

e”0T  *  t@-t.  (l3) 


This  identifies  an  exponentially  ssa.ll  neighborhood  of  t  within  which  the 
expansion  can  be  written  in  the  fora 

1/T  -  1/Tf  =  -x+2e“1itn(Tf/T)  +  . .. ;  (19) 

so  that,  to  leading  order, 

T  *  Tf/(l-xTf).  (20) 

Starting  at  the  value  for  x  =  0,  the  temperature  increases  without  bound 

as  t  increases  to  1/Tf.  The  unboundedness  is  a  consequence  of  our  failure 
to  account  for  reactant  depletion;  if  that  is  remedied,  T  increases  towards 
a  burnt  value  (see  equation  (25)),  entering  the  so-called  relaxation 

phase  when  it  is  0(9-1')  away.  The  relaxation  phase  lasts  an  exponentially 
short  time  also. 

These  features  are  shown  in  figure  3  and,  so  long  as  T  -  T.  *  0 ( 1 ) , 

D 

the  problem  without  depletion  provides  a  qualitatively  accurate  description 

% 

of  them.  In  particular,  the  fast  time  x  is  still  relevant. 

The  results  can  also  be  obtained  directly,  without  recourse  to  the 
exact  solution.  Thus,  if  the  expansion  (5a)  is  substituted  into  the 

* 

problem  (11),  we  find 

d$/dt  =  e*.  a  0  for  t  ■  0,  (21) 

with  solution  (IT);  and  introducing  the  fast  time  x  into  equation  (lla) 
yields 


with  solution  (19) • 


3(l/Tf-l/T-x) 

dT/dx  =  T*e 


(22) 


!+.  Inhoso geneous  Explosion 

We  now  combine  unsteadiness  with  spatial  inhomogeneity  by  considering 
the  slab  problem 

225  '9/T-8/T 

3T/3t  -3  T/3.X  *  (T^5/Tfe)(Tb-T)e  r  ,  (23) 


T  »  Tf  for  t  =  0  and  at  x  *  ±1.  (2*0 

To  account  for  reactant  depletion,  V  has  been  replaced  by  t7Y/Yf,  where  Y 
is  the  mass  fraction  (i.e.  concentration)  of  the  reactant  and  Yf  its  initial 
value.  The  Shvab-Zoldovich  relation 

T  +  Y  =  Tf  +  Yf  =  Tb  !25) 

(cf.  section  2.2)  is  then  used  to  eliminate  Y,  thereby  ensuring  conse.  >n 

of  total  enthalpy  (the  sum  of  thermal  enthalpy,  represented  by  T,  and 

chemical  enthalpy,  represented  by  Y). 

* 

During  the  induction  phase,  the  0(9-1)  departures  of  T  from  Tf 
expressed  by  the  expansion  (5)  satisfy 

3$/3t-3%/3x^*  <5e^,  ♦  =  0  for  t  *  0  and  x  *  ±1.  (26) 

Reactant  depletion  plays  no  role  during  this  initial  evolution  of’ the 

temperature.  When  the  system  is  subcritical,  i.e.  for  6  <  5c,  the 

perturbation  $  tends  to  the  steady  state  (8)  with  the  smaller  ,  as 

nr 

t  -►  *.  But,  for  6  >  $c,  the  absence  of  a  steady-state  solution  implies 
that  f  will  increase  without  limit,  and  indeed  thermal  runaway  is  found 
(numerically)  to  occur  after  a  finite  time. 

Further  progress  depends  on  a  description  of  this  runaway,  which  by 
symmetry  must  take  place  in  the  neighborhood  of  x  *  0.  Since  the  spatial 
derivatives  must  play  a  role,  they  have  to  be  increasingly  large  in  order 
to  be  comparable  to  the  ever-increasing  time  derivative.  It  follows  that 


-1.3- 


t he  region  in  which  runaway  occurs,  called  a  hot-spot,  must  continually 
shrink-,  this  self-focusing  is  an  essential  feature  of  the  process. 

The  appropriate  variables  for  the  runaway  are  t  and 

n  =  x/(te-t)J  (27) 

(see  figure  U),  where  t  (<5)  is  the  runaway  time,  to  be  determined  numerically. 
The  use  of  n  is  suggested  by  the  fora  of  equation  (26)  and  by  the  focusing 
discussed  above,  though  Kassoy's  achievement  in  identifying  it  is  in  no  way 
diminished  by  such  a  posteriori  observations.  Now  we  have 

3<j>/3t  -  (te-t)_1[»2/3n2-jn3!)>/3n]  =  e*  (28) 


and  we  seek  an  asymptotic  expansion 

P  =  -ln(t  -t)+t|>(n)+*  •  •  as  ttt  with  n  fixed, 

Q  Q 


finding 


(29) 


ip”  -  jniji’  +  e^  *  1  with  tp'(,0)  *  0  (30) 

(as  a  symmetry  condition).  Another  boundary  condition  is  needed  to  complete 
the  problem  for  ^ »  and  this  comes  from  matching  with  the  solution  outside 
the  shrinking  hot  spot,  i.e.  as  n  -*■  «  with  x  fixed.  Thus,  since  \p 
cannot  become  exponentially  large,  it  has  the  asymptotic  form 


ip  =  -2inn  A+...  as  n  -*■  +•,  (31) 

corresponding  to 

P  *  -2£nx-inS+A+. . .  as  t  t  t  with  x  fixed.  (32) 

Numerical  solutions  of  the  supercritical  problem  (26)  exhibit  the  behavior  (32), 
and  hence  determine  the  constant  A(5)  that  is  needed  for  the  second  boundary 
condition 

lim  (\|<+2inn)  *  A(5) . 


(33) 


The  problem  (30),  (33)  uniquely  3eterr.ir.es  the  function  j, ,  vhich  can  be 
readily  calculated  numerically. 

Following  the  initiation  of  thermal  runaway,  the  temperature  rises  at 
an  increasing  (exponential )  rate,  so  that  self-focusing  continues.  The 
variable  n  still  plays  a  role  during  this  process  but  now  the  appropriate 
time  variable  is  the  fast  one  (18);  in  terms  of  n,t  equation  (23)  becomes 

3T/3r+  3(jn3T/3n-32T/3n2)  =  (T2/Yf ) (Tb-T)exp(9/Tf-8/T-8T ) .  (2k) 

Note  that  x  =  6  5ne  with  n  =  0(1)  provides  a  measure  of  rapid 

focusing. 

The  homogeneous  result  (20)  is,  to  leading  order,  the  solution  here 
also;  otherwise  the  reaction  plays  no  role.  The  perturbation  is  found  to  be 

Tx  =  T^{ii(n)-ln[(l-Tft)(Yf-TfTbT)/Yf]}/(l-TfT)2,  (35) 

where  satisfies  the  equation  (30a).  It  must,  in  fact,  be  identical  to 
the  runaway  function  Just  constructed,  since  otherwise  there  would  be  a 
mismatch  as  t  -*■  0. 

The  hot-spot  evolves  so  rapidly  that  the  temperature  outside  has  no 
time  to  change,  i.e. 

T  =  Tf  +  8"1T2(-2tnx-ln6  +  A)+...  as  x+0,  (36) 

and  this  does  not  match  the  hot-spot  expansion,  even  to  leading  order.  The 
reason  is  clear:  the  expansion  (36)  breaks  down  at  points  inside  the  hot¬ 
spot  when  it  is  thickest  (t  small),  and  such  points  can  be  well  outside 
once  the  focusing  is  under  way  (t  moderate).  The  shrinking  hot  spot  leaves 
behind  an  intermediate  structure,  which  turns  out  to  be  stationary  (i.e. 
independent  of  t ) ;  to  leading  order  it  is  described  by 

T  -  Tf/(1-Tfx)  with  x  -  -20"1ln  x.  (37) 


-1.10- 


?or  the  homogeneous  problem  of  section  2,  vith  r.o  reactant  depletion, 
the  temperature  increases  indefinitely  as  t  -*  1/T^.i  but,  with  depletion, 
figure  3  shows  that  T  is  limited  by  tne  value  T^,  which  is  approached 
within  0(6  as  t  ■+•  (an  earlier  time).  The  latter  is  true  here 

also,  but  there  can  no  longer  be  the  single  relaxation  phase  shown  in 
figure  3  as  being  described  on  the  scale 


T  =  (t-tg)/£  with  £  *  (Yfe/Tpe 


-9Y./T,T. 


(38) 


since  now  the  temperature  is  close  to  only  in  the  vicinity  of  x  *  0. 

Instead,  there  is  a  transition  phase  at  the  hot-spot,  described  in  terms  of 
r  and  the  spatial  variable 


X  =  5sx/eJ , 


(39) 


followed  by  propagation  and,  finally,  relaxation  phases.  Figure  5  is  the 
result  of  numerically  integrating  the  limit  problem  for  the  transition  phase, 
during  which  the  focusing  of  the  hot-spot  is  opposed  by  reactant  depletion, 
thereby  forming  an  incipient  deflagration  wave  (section  2.U).  Once  formed 

the  wave  propagates  rapidly  through  the  right  side  of  the  slab  burning  up  the 

* 

reactant,  after  which  the  relaxation  phase  takes  over.  (Of  course,  similar 
remarks  apply  to  the  left  side  of  the  slab  too.)  Buckmaster  &  Ludford  (1982, 
p.  236)  give  the  details,  following  rCapila. 

5.  Ignition  by  External  Agencies 

So  far  we  have  confined  the  discussion  to  ignition  due  to  self-heating, 
but  it  can  also  be  caused  by  an  external  agency.  As  an  example,  suppose 
the  half-space  x  <  0  is  filled  with  a  combustible  material  subject  to  a 
prescribed  heat  flux  at  its  surface.  The  mathematical  problem  is 

3T/3t  -  32T/3x2  »  t?e"6y/T,j 


(2*0) 
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3i/3x  =  .^  >  3  at  x  3  0,  T  *  T^,  for  t  *  0  and  as  x  •*  -•.  (-*1) 
Here  the  constant  is  the  dimensionless  heat  flux  and  we  shall  suppose  that 

9/T 


V  3  «  r  '^th  Tr  >  Tf  ( i*2 ) 

The  carameter  T  characterizes  the  reactivity  of  the  material:  for 
*  r 

temperatures  below  T  ,  the  reaction  is  negligibly  weak  in  the  limit 

9  -*■  ®.  A  pre-exponential  factor  (even  depending  on  9)  can  be  given  to  V, 

but  this  is  equivalent  to  changing  Tr  slightly,  provided  nothing  is  added 

to  the  exponential  growth  of  V  with  9. 

During  an  initial  (finite)  time  interval,  the  material  is  colder  than 

the  reactivity  temperature  so  that  the  reaction  is  frozen  (i.e.  exponentially 

weak)  and  the  heat  equation  governs.  Because  of  the  heat  flux  T*  the 

s 

temperature  rises,  its  maximum  value  occuring  at  the  surface  x  =  0. 

Ignition  occurs  when  the  surface  temperature  reaches  T  .  The  subsequent 

r 

process  of  thermal  runaway ,  hot-spot  development  and  deflagration-wave 
formation  has  been  discussed  by  Kapila.  Here  we  shall  mention  only  the 
mathematical  problem  involved  in  the  thermal  runaway. 

The  rise  in  temperature  is  much  more  rapid  than  that  for  spontaneous 
combustion,  extending  only  over  a  period  0(9  ^).  If  conduction  is  to 
rival  the  temporal  changes,  spatial  gradients  must  therefore  be  0(9^)  so 
that,  at  corresponding  distances  from  the  surface,  the  temperature  has 

i 

dropped  by  0(9  !)  and  there  is  no  reaction,  (in  the  limit  9  -*•  «,  reaction 
only  occurs  at  temperatures  0(9  ^)  away  from  Tr.)  The  thermal  runaway 
is  therefore  governed  by  the  heat  equation 


3T2/3t  -  32T2/3i2  »  0 


(U3) 


and  its  cause  lies  in  the  nonlinear  boundary  conditions 
t+T2 

jTg/Sx  a  e  for  x*0»  as  X  and  as  t  *  (I*1*) 

Here  x  and  t  are  appropriately  transformed  space  and  time  variables  while 


T2  represents  the  second  perturbation  in  an  expansion  of  the  temperature. 

The  problem  (-2,43)  should  be  compared  with  that  for  spontaneous 
combustion,  given  by  equations  (26).  Clearly  they  have  entirely  differently 
forms . 

6.  Ignition  by  an  Externally  Generated  Hot-Spot 

Suppose  some  portion  of  an  infinite  combustible  material  is  burnt  very 
rapidly  so  that  there  is  a  local  rise  in  temperature  and  depletion  of  reactant. 
That  is,  we  create  a  hot  spot  (using  a  spark,  for  example)  somewhat  like  the 
one  that  develops  in  auto-ignition.  The  hot-spot  can  have  one  of  two  fates: 
either  it  decays  by  diffusion,  so  that  after  a  certain  time  the  temperature 
is  essentially  uniform  once  again  (and  constant  until  a  homogeneous  explosion 
occurs  because  of  self-heating) ;  or  it  acts  as  a  ignition  source,  producing 
a  deflagration  wave  that  sweeps  across  and  consumes  the  fresh  material.  If 
the  hot-spot  is  very  small,  temperature  gradients  will  be  very  large  and  the 
resultant  cooling  will  eliminate  it.  On  the  other  hand,  the  results  for 
auto-ignition  suggest  that  a  sufficiently  large  hot-spot  will  ignite  the 
material.  In  general,  the  initial -value  problem  that  must  be  solved  to  determine 
the  fate  of  a  particular  hot-spot  is  difficult.  We  shall  therefore  consider 
a  very  special  configuration  from  which  plausible  conclusions  can  be  drawn 
quite  easily. 

Consider  the  spherically  symmetric  form 


II  .  i_  -L  (r2  3T}  ,  v  (V*e/T 

3t  2  3r  Vr  3r '  v  '  Y_  e 
r  r 


(U5) 


of  equation  (23).  This  has  an  exact,  stationary  solution 

for  r  $  r„ 


h 

T-+Y  jr_/r 


(1*6) 
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in  the  limit  9  -*■  =°,  provided 


3/2  9/2T  2 

r#  *  Y^/2  9e  *f/20  T2. 


Thus,  the  reaction  term  is  asymptotically  zero  on  either  side  of  the  flame  at 
r  *  r#,  where  the  temperature  gradient  takes  a  Jump 


3T/3r j 


_ _  2  ,  -9/2T 

=  -  /2tf/Yf  T2  9"1e  b 


(cf.  the  deflagration-wave  solution  in  section  2.4). 

This  combustion  field  is  now  subjected  to  spherically  symmetric 

disturbances.  A  straightforward  stability  analysis  shows  that  the  perturba- 

t/Yf 

tion  of  the  flame  radius  grows  like  e  ,  which  corresponds  to  a  collapsing 
or  growing  hot  spot  depending  on  whether  the  flame  is  displaced  inwards  or 
outwards  initially.  The  result  suggests  that  the  radius  (46)  is  critical: 
larger  hot-spots  will  grow  and  smaller  ones  will  collapse. 
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Figure  Captions 


1.1  Heat-generation  and  heat-loss  curves. 

1.2  Steady-state  response  for  slab  with  surface  maintained  at  Initial 
uniform  temperature,  as  determined  by  equation  (9) . 

1.3  Temperature  history  for  homogeneous  explosion  with  reactant  depletion. 

1.4  Parabolas  n  ■  const. 

1.5  Temperature  profiles  during  transition  phase.  (Courtesy  A.K.  Kapila.) 
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